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In this work, we analyze the electroweak loop corrections to the Neutralino dark matter relic 
density in the framework of effective coupling. In the first part, we comment on the generic 
features of the corrections and quantitative changes to the predicted relic density. We analyze 
the correlation between the characteristics of effective couplings to the nature of neutralino. 
Effective couplings, however, absorb only the most dominant one loop corrections and are 
not an exact calculation. In the second part, we assess the validity of effective couplings by 
comparing them to the full one loop calculations in various regions of parameter space. 



1 Introduction 

Within the framework of ACDM, the existence of dark matter is well established. Although, 
the exact nature of this omnipresent dark matter is yet to be identified, its relic density is well 
measured and measurements from Planck satellite will reduce the uncertainly to 1%. This is good 
news because the abundance of dark matter not only depends on the astrophysical parameters 
like the Hubble constant and the velocity distribution, but also on the particle physics properties 
like it's mass and the cross-section. If dark matter is cold in nature then it's relic density Q, h 2 
is inversely proportional to the thermal average of the annihilation cross-section {av). Any 
changes in the theoretically predicted cross-section will thus change the relic density. In the 
light of upcoming results from Planck, a theoretically precise calculation of the relic density 
plays an important role in determining the favored regions of theory parameter space. 

Various theories beyond the Standard Model of particle physics (SM) seek an explanation for 
dark matter. Supersymmetry is one of the most elegant extensions of the SM not only because 
it provides us with possible dark matter candidate but also because it gives us an explanation of 
the hierarchy problem in the higgs sector Different ways of SUSY breaking at a high scale can 
lead to different dark matter candidates, the lightest neutralino or gravitino for example. Under 
the assumption of R-parity conserving SUSY, neutralino is the most widely studied dark matter 
candidate. However, the annihilation cross-section of neutralino is often taken at tree-level and 
changes due to loop-corrections can lead to a significant change in the predicted relic density. 
A full one loop analysis is numerically exhaustive due to large number of diagrams involved in 
annihilation cross-section while computing the relic density. 

In this work, we assess the impact of electroweak loop-corrections to neutralino dark matter 
via effective couplings and compare the validity of this approach to the full one loop analysis. 



2 Effective couplings for neutralino 



We implement effective couplings for the neutralino-fermion-sfermion vertex, as first discussed 
by Guasch et. al^ Analogous to the oblique corrections in the SM, they are flavor independent, 
finite corrections to the neutralino-fermion-sfermion vertex. 

The neutralino is a combination of supersymmetric counterparts of the gauge and Higgs 
boson, namely the gauginos and the higgsinos. After the SU{2)l x U(1)y symmetry breaking, 
the bino, wino and higgsino mix to form four neutralinos, the lightest one being of interest to 
us as a dark matter candidate. At tree level, the neutralino-fermion-sfermion couplings can be 
written as: 
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In the above equation, N denotes the neutralino mass diagonalization matrix, g is the SU (2) 
gauge coupling, 6\y the weak mixing angle, T3 and Qf are the third component of the weak 
isospin and the electric charge of fermion /, respectively, Pl,r = (1 ± 7s)/2 are the chiral 
projectors, mf and My/ are the masses of fermion / and of the W boson, respectively, and 
tan/3 is the ratio of vacuum expectation values of the two neutral Higgs fields required in the 
MSSM. Finally, the function //(/3) appearing in eq.([T]) is sin/3 for up-type quarks and sneutrinos, 
and cos/3 for down-type quarks and charged leptons; similarly the index hj = 4(3) for up-type 
quarks (down type quarks and charged leptons). 

At one loop level eq.([l]) receives counter-terms and one can extract a flavor independent set 
of corrections to the neutralino mixing matrix as follows: 
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Here, 5g, St\y, 5M^r, 5 cos /3, S sin/3 represent the counter-terms of the quantities introduced in 
eq.(pQ) and SZ^f represent the neutralino wave-function renormalization constants. 

Note that eqs.([2]) include only two point correlation functions and hence do not include any 
vertex or box corrections. Also, it should be remembered that these corrections modify mixing 
matrix N appearing only in the interaction lagrangian. 

Eqs.© are finite only for the fermion-sfermion loops. At this point a couple of remarks are in 
order. First, the above expressions are process independent in the sense that they do not depend 
on the details of external state, second, they were originally introduced to absorb the one loop 
corrections in the neutralino-fermion-sfermion vertex, however, they can also be used to (partly) 
absorb all the vertices involving the gaugino component of the neutralino, except for the vertices 
involving gauge bosons. The most interesting feature of the above set of effective couplings is 
the non-decoupling behavior of heavy sparticles. Due to broken super symmetry, otherwise equal 
sfermion-fermion-gaugino coupling g runs apart from the corresponding fermion-fermion-gauge 



boson gauge coupling g below the sfermion mass scale. At one loop level, this difference is given 
by: 
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Thus, correction to g logarithmically increases with the sfermion mass scale. 
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3 Results 



The implementation of these effective couplings is done in the on-shell scheme. We use a variant 
of the on shell scheme for the neutralino sector, as introduced in Chatterjee et.alJ^We define the 
electromagnetic coupling constant, W and Z-boson masses as on shell quantities. To calculate 
the relic density, we modify the couplings in micromegas with the above expressions. Taking 
the input at electroweak scale, we scan over the bino parameter M\. We set Ma = 500 GeV, 
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Figure 1: Correction to the bino-like component of the neutralino mixing matrix (left) and the corresponding 

percentage change to the relic density (right). 

we plot the correction to the bino component of the lightest neutralino AN±\ (left) and the 
percent change in the relic density Af2 (right). Since M2 = 400 GeV, for M\ > 400 GeV, the 
lightest neutralino is not bino-like anymore. As a result, the correction AiVn decreases. While 
the neutralino is bino-like, the changes in the relic density are characterized by the corrections 
captured by AiVxi and opening of various annihilation channels, this is reflected by various peaks 
appearing in the plot. In some of the annihilation channels only part of the radiative corrections 
are captured (for example, at M\ ~ 84, 94, 106 GeV, the W + W~ , ZZ, Zh channels mediated by 
a Higgs exchange open and the radiative corrections only to the gaugino-higgisino-neutral biggs 
vertices are captured.). 

4 Comparison to full one loop analysis 

This being an effective approach, will fail at some point and it is important to know in which 
regions of parameter space such couplings are applicable. To this extent we compare the results 
obtained for cr(xiXi ~* ^ + using effective couplings to those of full one loop results. ^ The 
full one loop results were obtained by using SLOOPS. 

In the case, we set M A = 1 TeV, tan/3 = 4, Mj = = 500 GeV, M fifl = = Mq l = 
800 GeV, A f = 0, M 2 = 500, M 1 = 600 GeV and fi = -600 GeV. In Fig. El left hand side 
shows the percent change in <t(xiXi ~ * ^ f° r the full one loop calculations, the effective 
coupling approach and the difference between the two. We scan over M\ changing the bino 
fraction of the neutralino. The full one loop analysis shows about 18 to 20% corrections of 
which about 15 to 17% corrections are absorbed by effective couplings remaining non-effective 
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Figure 2: Comparison of the full one loop analysis with the effective coupling results, on left, a scan against the 
bino-like component of neutralino Mi and on right against the higgsinolike component fi. 

part staying constant at about 2 to 3%. This is because the neutralino stays bino-like and 
annihilates dominantly to fermion final states via a sfermion exchange, thus the box and vertex 
corrections do not contribute significantly. The figure on right shows the percent change in the 
same cross-section as a function of the higgsino component of the neutralino [i. In this case, the 
neutralino annihilates via a Z boson exchange and the box and vertex corrections are significant. 
Effective couplings fail to reproduce the full one loop correction, which ranges from to — 40%. 
The green line, illustrating the contribution of the boxes is the most dominant contribution, 
which the effective couplings do not take into account. We tried to implement a new set of 
effective couplings, displayed in yellow line, however even they fail to reproduce the genuine box 
diagrams. 

5 Conclusions 

Effective couplings as used in this work can be an easy and convincing way to include dominant 
process independent one loop electroweak corrections to neutralino dark matter relic density. 
These are constructed using two point correlation functions and can be used to absorb corrections 
to all the gaugino component of the neutralino. When applied to a dominant bino-like neutralino, 
the changes to the relic density can be as large as 4%. They are most effective when the neutralino 
is bino-like, annihilating into a pair of fermions via sfermion exchange when the contribution 
from the box and the vertex corrections are negligible. When the neutralino changes it's nature 
to wino or higgsinolike however, the effective coupling approach fails to reproduce full one loop 
results. 

6 Acknowledgements 

I thank the Moriond organizers for financial support to participate in this conference. 
References 

1. M. Drees, R. Godbole and P. Roy, Hackensack, USA: World Scientific (2004) 555 p 

2. J. Guasch, W. Hollik and J. Sola, JHEP 0210, 040 (2002) [hep-ph/0207364] . 

3. A. Chatterjee, M. Drees, S. Kulkarni and Q. Xu. larXiv:1107.5218l [hep-ph]. 

4. F. Boudjema, G. Drieu La Rochelle and S. Kulkarni, Phys. Rev. D 84 (2011) 116001 
[arXiv:1108.429T1 [hep-ph]] . 



"Note that we set a em at a ern (q 2 — 0) in this case as opposed to results of Fig.[T]where a em is at a em (q 2 = M§) 



